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Abstract 

The cosmological evolution of a quintessence-like scalar field coupled to matter and 
gauge fields leads to effective modifications of the coupling constants and particle masses 



c/3 ■ over time. We analyze a class of models where the scalar field potential V{(f)) and the cou- 

plings to matter -8(0) admit common extremum in 0, as in the Damour-Polyakov ansatz. 



We find that even for the simplest choices of potentials and 5(0), the observational con- 
straints on Aa/a coming from quasar absorption spectra, the Oklo phenomenon and Big 
Bang nucleosynthesis provide complementary constraints on the parameters of the model. 
We show the evolutionary history of these models in some detail and describe the effects of 
a varying mass for dark matter. 



1 Introduction 

The existence of dark energy supported by a number of cosmological observations remains 
one of the most outstanding puzzles in modern cosmology [1]. The cosmological constant 
and/or a quintessence field are the most commonly proposed candidates for dark energy. 
The possibility that a scalar field at early cosmological times follows an attractor-type 
solution [2, 3] and tracks the evolution of the visible matter-energy density while dominating 
the energy density in the Universe at late cosmological stages has been a subject of lively 
debates in the cosmology literature [4] - [20] . There are certain hopes that quintessence-like 
models may help alleviate the severe fine-tuning associated with the cosmological constant 
problem. This by itself represents a legitimate field of research and has triggered various 
model building efforts. 

The most fundamental difference between quintessence and the cosmological constant 
is that the former represents a new very light degree of freedom, with a "wavelength" 
comparable to the size of the Hubble horizon. As a consequence, the dynamical scalar 
field gives rise to a time-dependent equation of state parameter, oj^ = p^/ p<f,. Therefore, 
the variation of u for dark energy with redshift is among the most important cosmological 
parameters to measure, and will be a large component of future efforts in cosmology. 

Are there additional experimental ways of checking for the existence or absence of dark 
energy in the form of quintessence? Apart from "geometric" tests using standard cosmolog- 
ical methods (CMB, supernovae, lensing, etc.), one could also hope to detect the possible 
interaction of the quintessence field with matter. This can occur either through the ob- 
servation of the effective change of the coupling constants and masses of particles over 
cosmological times, or via detecting additional components to the gravitational interaction 
due to exchange by the scalar. A positive detection of the cosmological change of coupling 
constants would be firm proof of the existence of new degree(s) of freedom with extremely 
long wavelength, thus providing a perfect candidate for quintessence. 

For this reason, the recently reported indications of a change in the effective electromag- 
netic coupling constant [21], Aa/a ~ —0.6 x 10~^ at redshifts z ~ 1, if indeed true, can be 
considered as an independent detection of an ultra-light degree of freedom. These results 
[21] are based on the comparison of the relativistic shifts of atomic energy levels in quasar 
absorption spectra. For some time this result remained unchallenged, yet recently there 
have been attempts to detect a variation in a using similar methods [22] that has so far led 
to null results. Also, alternative (a-unrelated) explanations of the results of Ref. [21] based 
on a z-dependent Mg isotope abundance were shown to be astrophysically viable [23]. This 
rather controversial situation becomes even more complicated if other constraints on Aa/a 
are taken into account [24]. The Oklo natural reactor and meteoritic abundances of rhe- 
nium provide stringent constraints on the change of the coupling constants that goes back 



to z ~ 0.1 — 0.4 [25, 26, 27]. However, a recent re-analysis of Oklo phenomenon suggested 
that the present data are consistent with the non-zero change of Aa/a = 4.5 x 10^* [28]. 

Despite the questionable status of the non-zero claim for Aa/a there have been a 
number of attempts [14] - [17], [29] - [35] to build simple models that could account for a 
possible O(10~^) relative shift in the fine structure constant at redshifts 2; ~ 1. It is widely 
accepted now that the minimal Bekenstein model [29] with a scalar field coupled only to 
the electromagnetic field in the initial Lagrangian, 

^- - -> {^) -^^-" -K^ ^ ''^. ^'ii^S^ ■■) "-"" • <'■'' 

cannot provide Aa/a larger than 10^^ — 10^^° level when the constraints on the nonuni- 
versal gravitational interaction mediated by this field are imposed. This is because the 
electromagnetic portion of the baryon energy density that drives the cosmological evolu- 
tion of (/9 is a very weak source. It has been suggested that the coupling to dark matter 
[30, 35, 33] and/or the self-interaction potential [34] of the scalar field drives its evolution. 
A number of tracker field and quintessence field models coupled to F^^ have been analyzed 
with the universal conclusion that a linear coupling to the electromagnetic energy density 
on the order of 10~^ ^ (f ^ 10~^ can easily explain the QSO-suggested change of the fine 
structure constant and satisfy the experimental limits on the universality of the gravita- 
tional interaction. It has proven to be a harder task to have a substantial change in the fine 
structure constant at 2; ~ 1 and remain consistent with the Oklo and meteoritic constraints. 
Only a few models (see e.g. Refs. [15], [17]) are known to pass these requirements. 

To this end, it will be useful to investigate a class of quintessence models where the 
potential V{ip) admits a minimum at ipo, and ip{z) approaches <^q when z ^ 0. Without 
loss of generality, we can put the extremum value to zero, and then ^(0) ~ A = I^aPc where 
Pc is the closure energy density today. In order to suppress the variation of a at late times, 
one could also assume that (f = 0, or in other words, require V{ip) and Bf{ip) to share the 
same extremum. We note that most of the quintessence models analyzed in the context 
of changing coupling constants to date, were coupled to the electromagnetic lagrangian 
linearly, while the coupling to other gauge and matter fields was neglected. Continuing 
along the same line, it would be reasonable to allow a multitude of couplings Bi{ip), and 
require Q = 0. This is exactly the proposal of the "least coupling" principle, introduced a 
decade ago by Damour and Nordtvedt [36] and Damour and Polyakov [37] in the context 
of the string dilaton with the primary goal to suppress the strength of the gravitational 
interaction mediated by ip. Here we supplement their idea by including the self- interact ion 
potential V{ip) with an extremum at the same value of ip as in Bi{ip). We note that this 
property, a common extremum in ^9, remains intact even in the presence of the radiative 
corrections, and therefore does not require any additional fine-tuning. 

In what follows, we analyze in detail the cosmology of coupled quintessence models with 



a common extremum in V{ip) and gauge and matter /gauge couplings Bi{ip). We note that 
only a coupled quintessence model derived from Brans-Dicke theory has been analyzed to 
date [38, 39]. We investigate the size of the possible variation of masses and couplings in 
these models and find that some generic choices of V{(f) and Bi{ip) can be consistent with 
all observational requirements. This paper is organized as follows. In the next section we 
introduce our model and display the necessary cosmological equations. In sections 3-5, we 
make specific choices for the potential V{ip) and matter /gauge couplings Bi{ip) and find the 
evolution of the dark energy and dark matter energy density over the redshift. We make 
predictions for the variation of the coupling constants and masses and outline the choice of 
parameters that satisfies all existing constraints. We reach our conclusions in section 6. 

2 Cosmological evolution of a scalar field in the pres- 
ence of couplings to matter 

We start by writing down the general equation for the interaction of a light scalar field (/? 
with matter, 

s^ = J d^xV^l—[dyd,(i) -R]- v{<p) 

_M^F»F«'^'^ + Efe^^. - B,{<p)m,ij,ijA (2.1) 

j ^ 

In this expression, is the scalar field of interest with its kinetic term normalized together 
with the gravitational interaction. The placement of the reduced Planck mass, M, in 
(2.1) renders dimensionless. Bpi^cp) represents the 0-dependence of the gauge couplings 
in Standard Model where the sum is over all three groups, ■j/'j represents Standard Model 
fermions that are coupled to via the functions Bj{(f)). Note that the 0-dependent rescaling 
of the matter fields and metric allows us to remove from the couplings to the kinetic terms 
and R. A similar rescaling in the gauge sector will lead to the appearance of a 0-dependence 
in the interaction terms between gauge fields and matter, i.e. precisely 0-dependent gauge 
couplings. Besides SM fermions, the sum over j includes other matter field {i.e. scalar 
Higgses, Majorana neutrinos etc.) as well as various interaction terms. Finally, any viable 
cosmological scenario requires cold dark matter, and we also include it in the sum with the 
separate coupling B£)m{<P)- Since couples to the trace T^ of dark matter, our results are 
independent of the nature of the dark matter particles (scalars or fermions). Furthermore, 
since dark matter is the dominant component of the matter energy-density, we can simply 
take 5„(0) ~ Bdm{,<P)- 



From the above action (2.1) we can derive the equation of motion for the scalar field, 0: 

where (ipjipj) stands for the number density of a j'-th fermion. Notice that at tree level, 
radiation does not contribute to the r.h.s. of Eq. (2.2) because its stress-energy tensor is 
traceless. Given a potential V{(f)), the perfect fiuid energy density and pressure contribu- 
tions due to are: 

p^ = ^M^' + V{(i^) (2.3) 

P^ = \m^''-V{<P)^uj^p^, (2.4) 

where the parameter uj^ is related to the equation of state (EOS) of the scalar field. Com- 
bining Eqs. (2.3) and (2.4), we obtain a useful relation between the potential V^(0), the 
energy density of the scalar field, and a;^, 

^ - 2V(0) p.5) 
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The energy density p<^ and pressure p^ of the scalar field contribute to the r.h.s. of Einstein's 
equations and yield the following Friedmann equation in a Robertson- Walker Universe, 



^ = 3M^ V^^ + ^'^ + M - ^M^^-' (2.6) 

where H = a/ a is the Hubble expansion rate, M^ = M^/Stt is a reduced Planck mass, 
Pr and pm is the energy density of radiation and matter, and pcr is the critical energy 
density. The energy density of matter and radiation contains ^-dependence, i.e. pm = 
J2j Bj{(f))mj{iljjiljj). Using these definitions, we can rewrite the scalar field equation (2.2), 

d) + 3H6 + ^-^— = — ^ — TT^^Pm- (2.7) 

^ ^ M2 90 M2 90 '^ ^ ' 

Conservation of the energy-momentum tensor gives us another equation, 

p^ + 3H{l + uj^)p^ = 'baT'P'^^' ^'^'^^ 

For cosmological studies that span a large range of redshifts z, it is convenient to intro- 
duce the variable x as the logarithm of the scale factor a, 

a; = lna = -ln(l + ;z) (2.9) 



where we choose the present scale factor a^^' = 1. With the use of the variable x, we can 
rewrite the relevant system of equations for dlnpi/dx in the following form, 

dlnpm _ „.^ ^ din Bm{(l)) d(p 

dlnpr 



dx 



3(1 + CJ,), (2.11^ 



dlnp^ ^ _ d\nB^{(f)) d(f)p^ 

dx d(f) dx p^ ' 

where Ur = 1/3 and c^m = should be used for radiation and matter respectively. 

The evolution of as a function of redshift can be found relatively easy if one is able 
to solve the cosmological equations and obtain u^ and Qi = Pi/pcr as a function of x. For 
example, the derivative of with respect to x can be obtained from Eq (2.3), 




-f] =3n^{l + uj^), (2.13) 



If u!^ is close to —1, the kinetic energy of the scalar field goes to zero which occurs when 
the scalar field is close to the minimum of the potential. 

The evolution of 0(a;) expressed in terms of uj^ and ^2^ takes an especially simple form 
when Bm{(f)) does not produce a significant change in the cosmological equations, and 
terms containing dBm{4>)/d(f) can be neglected. In this case, the matter and radiation 
energy densities can be easily related, pr = {aeq/a)pm, where aeq is the scale factor when 
the radiation and matter densities are equal, and a^q ^ 1. This allows us to rewrite Eq. 
(2.6) for the evolution of the scalar field energy density in a useful form. 



dx 



V (7, + (7, _ _ / 



Ct' I lA/flQ 



(2.14) 



The dependence of critical density on the scale factor can be found explicitly, which allows 
us to express V{(f)) from (2.5) as 

V{<P) = 1(1 - uj,)Q,p,, = ^(Mi7(°))2(l -.;,)fi,i^ (v)'^^- ^^-^^^ 

In this expression, H^'^^ and I — Vli = f2^^ are the Hubble expansion rate and the matter 
density today. This equation can be used further to obtain 0, i.e. by taking the logarithm 
of both sides of (2.15) for a simple exponential potential. The derivative of (2.15) with 
respect to x gives us another useful equation, 

Tx = '^' + "^) + ^^^- ^2.16) 



When the change in Bm{4>) is not small and cannot be neglected, all cosmological equa- 
tions become considerably more complicated. The scaling of the matter energy density 
differs from usual a~^ behavior because of the changing mass, due to Bfn{4>), while the 
scaling of radiation energy density remains unchanged. 



Pmix) = Pl^'^ t" !!!!' ' ^^^ P.(:r)=p(°)a-^ (2.17) 

5„ 0(0) 



Pr ^TTTTT^P™ K^-^^) 

a bJMxY 



(2.19) 



With the use of these relations, we generalize Eqs. (2.14) and (2.16), 



dx 






a + at^^\ ^a + aij^ d(p dx 



d In Bm 



dln{l-uj^) dlnVdc/) (I - il^) ^ a ^^\nBrnd(j) 



dx (90 dx n^ ^a + afq ^^ dx 

where we have introduced an auxiliary function, a'f^{(f)) = aeqBrn{(f)ixeq))Bj:^^ 



(X 



3 Cosmological evolution of the scalar field 

3.1 Relation between V{(j)) and Bi{(f)) 

Our main assumption in this work is that all functions Bi{(f)) and V{(f)) admit a common 
extremum, which is an obvious generalization of Damour-Nordtvedt and Damour-Polyakov 
constructions [36, 37]. Near this extremum, all functions admit an expansion 

5,(0) = 1 + 1^.0' + ■••; v{<p) = \/o(l + ^A02 + ...), (3.1) 

where C,i and A are dimensionless parameters, while Vq is of the order of the dark energy 
density today. 

In order to be able to start our analysis, we must specify the potential V{(f)) and functions 
Bi{(f)). Without any serious guidance from the underlying theory, their choice is completely 
ad hoc, and the only constraints are those coming from observations. A large number of 
quintessence potentials V{(j)) has been analyzed in the past, but the parameter space of 
functions Bm{4>) and Bp{(f)) remains relatively unexplored. To this end, we propose an 
ansatz that would relate Bi{(f)) and V {</)). Since we are going to consider only positive- 
definite functions, we choose a two-parameter relation, 

B,i^) ^ (h±ml}i^\"\ where6. + l>0 (3^2) 



This relation is sufficiently generic, and varying the two dimensionless parameters rii and 
hi for each function Bi{(f)) allows us to cover a wide range of possibilities. The expansion of 
(3.2) near the extremum allows us to relate C,i and A, 

Xrij 

i. = YTt- (^^^' 

If the quintessence field changes by A0 ~ 0(1) between redshifts 2; ~ 1 and z = 0, the 
constraints on the variation of masses and couplings would normally require small C,i which 
may be due to small Ui and/or a large value for 6j. In a typical tracker regime, V{(f)) changes 
by many orders of magnitude, and at early cosmological epochs, the value of V{(f)) can be 
much larger than the normalization scale 6j. In this case, Bi{(f)) ~ (1^(0))". 

In the next two sections, we analyze two types of potentials, which in the asymptotic 
regime of tracking behavior, scale as exp(0^) and exp(0). We ffist assume that the couplings 
of to all matter fields are rather small and have no backreaction on the cosmological 
evolution of 0. We then allow for a substantial coupling Bm{4>) to dark matter and show 
the impact of that on the dynamics of the quintessence field. By comparing the evolution of 
-Bir(0) for electromagnetic field to observational constraints, we can determine the allowed 
range of parameters for np and bp- 

3.2 Quintessence with exp(A0^/2) and cosh(A0) type potentials 

We next consider two sample potentials which can combine the properties of a tracking 
solution with a minimum that provides the dark energy at late times: 

case A: \/(0) = \/oexp(-^) (3.4) 

case B: V{(f)) = \/ocosh(A0) (3.5) 

Both potentials have the same expansion near the minimum (3.1) and are sufficiently easy 
to analyze. For now, we set Bp = Bm = 1. Closely related potentials have been expansively 
analyzed in the past by several groups [7, 9]. 

Case A. 

While the evolution equation for i7^ is unchanged from Eq. (2.14), the evolution equa- 
tion for cj0 can be rewritten as 

d\n{l — ujfj,) 



3(l + cc;^)±A0y'3fi^(l + cc;^), (3.6) 



dx 

where the choice of plus(minus) depends on whether the evolution of starts from nega- 
tive(positive) values. The dependence of can be obtained directly from (2.15), 

"3(ijWM)2(l - nf] 



\^^ = -4x + ln(a + a,,) + In (^^^^) + In 



2V. 



(3.7) 



All terms in Eq. (3.7) other than the (— 4x)-term do not change significantly over the 
history of Universe. This is similar to the evolution of a scalar field in the tracker regime 
for an exponential potential, where lnV{(j)) ~ const — 3(1 + uji)x [11]. 

An advantage of a tracker-type potential is that the cosmological evolution of the scalar 
field during the observationally relevant epoch, i.e. z < 10^° is insensitive to the choice 
of initial conditions for that can be specified deep inside the radiation-dominated epoch. 
In order to exhibit tracking behavior, the potential must satisfy two conditions [4]: 



V"V 



rfln(r-l) 



dx 



< 1. (3.8) 



The simple exponential potential provides a tracker solution with F = 1. 
For the potential (3.4), we find 

1 



A02' 
which is larger than 1 for any positive A. For the second condition we obtain 

dln(F-l) _ 2^3Q^{l + uj^) 



(3.9) 



dx 



(3.10) 



which is much smaller than one as long as is far from the minimum, i.e. at all cosmological 
epochs except for late times. 

To be compatible with observational data, quintessence models must satisfy the fol- 
lowing conditions: 1) The energy density of quintessence must be subdominant during 
Big Bang nucleosynthesis (BBN) [3], nf^^\x ~ -23) < 0.2 at T ~ 1 MeV; 2) the 
energy density at present must be compatible with the preferred range for dark energy, 
0.60 < f^^ (x = 0) < 0.85 [40]. The evolution of energy density components, ilr,^m, and 
Q^ is shown in Fig. la along with the dark energy equation of state. In the very early 
universe, the EOS parameter has an oscillatory behavior as does the density parameter Q^p 
which quickly converges to an attractor solution washing away any memory of initial con- 
ditions. As one can see, the constraint of Q^ < 0.2 from nucleosynthesis is easily satisfied 
in this model. We also note that during the radiation dominated era, u;^ ~ cu^ = 1/3 as 
in other tracker solutions. The transition to a matter dominated Universe occurs around 
X = —8.6 {z ~ 5400). During the matter dominated stage, u^ begins to drop sharply 
towards its final value of u^ = —1. At x = —0.5 {z ~ 0.6), the Universe begins to be dom- 
inated by the quintessence field. In Fig. lb, we show the present value of tu^ as a function 
of the coupling. A, for three different values of the present value of Q^. For large values of 
A, the potential is relatively steep and the field quickly evolves towards its minimum where 
\/ ~ Vo, which must be tuned close to the present cosmological constant value. In this case, 
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Figure 1: The exp(A0^/2) potential, a) The cosmological evolution of the equation of state 
parameter, tu^, and the energy density parameters, fii, of each component for A = 5. The 
evolution of the parameters is similar for other choices of A. b) The dark energy equation 
of state parameter uj^ as a function of A for an acceptable range of the dark energy density 
(0.6 < l^Jf^ < 0.85). 

u^ ^ —1. In contrast, for small values of A, the potential is relatively shallow, and is still 
evolving. As one can see, applying the constraint [41], lu^ < —0.6 gives the limit A > 1.5 
(for the case of fi^ = 0.73). When the present value of is not at the minimum (0 = 0), 
the constant, Vq must be adjusted to obtain the desired value of fi^ today. 



Case B. 

As one might expect, the potential (3.5) gives a very similar result as one would obtain 
for a simple exponential potential (V^(0) oc exp(— A0)) at early times. Therefore, for A0 > 1, 
i.e. far from the minimum, the scalar field is simply 



A0 ^ —Ax + ln(a + aeq) + In 



1-a 



+ ln 



3(i7WM)2(l - nf) 
2^ 



For this potential, the tracking parameter F is. 



1 + 



sinh (A0) 



(3.11) 



(3.12) 



Clearly, at late times, as approaches 0, the tracking condition will be violated which 
opens up the possibility that the Universe will be dominated by the scalar field, in contrast 
to the simple exponential potential. A comparison of the late time behavior between the 
exponential and cosh type potentials is shown in Fig. 2. 
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Figure 2: The evolution of Q^p and u^ for the simple exponential potential (solid) and the 
cosh potential (dotted). The main difference between two potentials appears at late times. 
We use A = 5 for both potentials. 



In Fig. 3a, we show the full evolution of Qr,^m, ^ip, and Ufp. In this case, matter 
domination occurs at x = —8.2 {z ~ 3500), and quintessence domination occurs at x = —0.3 
[z ~ 0.35). The cosh-type potential turns out to be more restricted by observational 
requirements than the potential (3.4). In Fig. 3b, we show the value of Q^{x = —23) 
corresponding to its BBN value as a function of the coupling A. We see that the constraint 
(^0 < 0.2 implies 4 < A. 



4 Time variation of the fine structure constant 



To study the cosmological evolution of the fine structure constant we use the following 
relation between Bp and a, 



Aa{z) _ a{z) - a(0) 5ir((/.(0)) 



1. 



a 



(4.1) 



a(0) BF{<f){z)) 

Note that the present value of the field 0(0) is close to zero. 

We assume that the backreaction of Bi{(f)) on the cosmological evolution of itself is 
small, and therefore can use the results of the previous section. All physical processes 
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Figure 3: The cosh(A0) potential, a) The cosmological evolution of the equation of state 
parameter, u^, and the energy density parameters, Qi, of each component for A = 5. The 
evolution of the parameters is similar for other choices of A. b) The quintessence density 
parameter Q^ evaluated at the time of BBN as a function of A. 

sensitive to changes in a can be separated into two broad categories: those that correspond 
to relatively low redshifts, z ~ 1 and smaller, and high-redshift phenomena such as BBN 
and the cosmic microwave background anisotropics. 

4.1 Late time evolution 

For the late time evolution of a we can use the expansion of Bf{4>) given in (3.1). In this 
case, the expression for A (a) takes the following simple form. 



Aa(z) 



a 



I {<P\0) - <p\z)) 



(4.2) 



where we dropped the subscript F in C,f to be concise. For the exponential of 0^ potential 
(3.4), this expression can be found in a more explicit form with the use of (3.7), 



Aa(x) 



a 



A 



)(0)^ 



-3a; + In 



(]^(l-fi7)(l-a;^) 



>(o) 



(0)n 



^r(i-^<^)(i-^r) 



(4.3) 



Using the result of the previous section, we can predict the evolution of a over redshift 
in terms of two parameters, ^ and A. We choose two characteristic values of ^, based on two 
QSO results. To be consistent with the non-zero result for Aa by Murphy et al. [21], we 
choose ^ in such a way that Aa/a = —5.4 x 10~^ at a redshift of 3. Another option that we 
explore is |Aa/a| < 6 x 10^^ at ^ = 1.5, which is motivated by the experimental accuracy 
of Chand et al. [22]. For definiteness, in the second case we choose Aa/a = —7 x 10^^. 
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The results for the choice l^aja = —5.4 x 10^^ are plotted in Fig. 4. The family of solid 
curves corresponds to different choices of the coupling A = 1, 4, and 5 in (3.4). It is easy to 
see that the evolution of a is approximately linear over a large range of redshifts for z > 1. 
This follows directly from (4.3) where only the —3a; term is changing at these redshifts. 
In panel b) of Fig. 4 we show a blow-up of panel a) for z < 0.5. This will facilitate the 
comparison between this model and the experimental constraints at low z. At late times 
(at low z) , the evolution of a slows down considerably. For comparison, we present similar 
results for the case of the linear coupling ( between and F^^.F'^'^, for the same choice of 
initial conditions. The corresponding curves are plotted in Fig. 4 by dashed-dotted lines. 
The main qualitative difference between these two choices, ^ or (^ is at late times. When 
C = and .^ 7^ 0, the fine structure constant evolves very slowly, whereas for ( ^ the 
evolution of a remains approximately linear. 



In Fig. 5, we show the analogous results for the evolution of a for the potential (3.5). 
Once again, we show results for three choices of A; in this case, we have used A = 4,4.5, 
and 5. Notice that in this case, a overshoots its present day value. This is due to the fact 
that with these choices of A, (p rolls past the minimum (at = 0) and begins an oscillatory 
behavior as it settles to the minimum of the potential. For the quadratic coupling, variations 
in a do not change sign thus causing Aa to reach the value 0, before becoming negative 
again. Damped oscillations of (and hence a) will continue to settle towards the minimum 
for which Aa = 0. 



We now turn to the question of how the sensitivity to Aa/a achieved in Refs. [21, 22] 
compares to other probes of Aa/a at smaller redshifts. In particular, we calculate Aa/a 
dX z = 0.14 relevant for the Oklo bound [25], the time average Aa/a between the redshifts 
and z = 0.45 which is constrained by meteoritic data [26, 27], and the rate of change 
in a today, a/ a [45]. In addition to the observables directly related to the fine structure 
constant, we also calculate the differential acceleration of two bodies towards a common 
attractor, limited by the precision tests of the universality of the gravitational force [35, 46]. 
To calculate Ag/g we use the calculation of Ref. [35], and conservatively assume that the 
coupling of the scalar field to nuclei is mediated only by the electromagnetic interactions. 
The calculations of Ref. [35] are done for the linear coupling (, but the generalization to 
the case of C = 0,.^ 7^ is straightforward, via the substitution ( -^ C,(f){z = 0) [37]. The 
results of these calculations are compiled in Table 1 . 

In Table 1, we display specific results for both forms of the potential, V {</)), considered; 
both with A = 5. We also show results for both sets of normalization [21, 22] and for 
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Figure 4: a. The evolution of Aa/a over the redshift range < 2; < 3 driven by the 
potential (3.4). Panels a) and b) use the common normalization Aa/a = —0.54 x 10^^ at 
z = ?). Figures c) and d) use the common normalization Aa/a = —0.06 x 10~^ at z = 1.5. 
The solid lines correspond to the choice C = ^-^d ^ 7^ 0, whereas the dashed-dotted lines 
allow C 7^ 0. 
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Table 1: Results for /\a/a for comparisons with the Oklo, meteoritic, present-day, and 
equivalence principle bounds. Values of iS^aja^ ^, and ( have been scaled by a factor of 
10^, values of a/a and Ag/g have been scaled by a factor of 10^^. 



Vo 


(^)3 


(^)l.5 


e 


C 


( a )o-14 


(Aa\ 


f (yr-^) 


A.9 
.9 


exp(A0V2) 


-5.4 


-3.4 





-7.6 


-0.24 


-0.42 


9.6 


-2.5 




-5.4 


-2.1 


22 





-0.011 


-0.05 


0.048 


-0.074 




-0.93 


-0.6 





-1.3 


-0.041 


-0.073 


1.7 


-0.076 




-1.4 


-0.6 


5.8 





-0.0028 


-0.013 


0.013 


-0.0053 


cosh(A0) 


-5.4 


-0.54 





-29 


-0.097 


0.20 


0.17 


-36 




-5.4 


-0.054 


310 





-0.0017 


-0.029 


0.0014 


-1.5 




-4.5 


-0.6 





-24 


-0.082 


0.17 


0.16 


-26 




-34 


-0.6 


1900 





-0.011 


-0.18 


0.012 


-59 



both a linear {( ^ 0) and quadratic {^ ^ 0) coupling of to F'^. Here, we will not 
consider the constraints on a when relations between all three gauge couplings are assumed 
[42, 43, 44]. The latter are typically a factor of 100 times stronger. The Oklo bound [25] 
is roughly Aa/a < 10^^ Therefore any result in Table 1 with a value of (— )o.i4 in excess 
of 0.1 is excluded. For the exp(A0^/2) potential, the linear coupling with the —0.54 x 10^^ 
normalization is excluded, while the other cases are all allowed. From Fig. 4, we see 
also that smaller values of A lead to larger variations in a at 2; = 0.14. For the cosh(A0) 
potential, all cases with A = 5 are allowed. However, as one can see from Fig. 5, even the 
slightly lower values of A shown lead to excessive variations in a for the linear coupling. 

A closer look at Table 1 reveals an interesting observation. A non-zero value for Aa/a 
at the level of ~ 0.5 x 10~^ ai z = 3 can be perfectly consistent with the bound \Aa/a{z = 
1.5)1 < 0.6 X 10~^ because of the highly non-linear evolution of the scalar field. Indeed, the 
quadratic coupling for the cosh- like potential shows that the ratio Aa{z = 1.5)/Aa{z = 3) 
can be as small as 0.01. This implies that in certain models the non-zero result of Ref. 
[21] for Aa/a that spans a rather large range of redshifts and the better-sensitivity zero- 
result measurement of Aa/a at one redshift z = 1.5 may not necessarily be contradictory. 
This would require, however, that the dominant source of the effect observed in [21] is 
due to their high redshift absorbers {z > 1.8) which also carry the largest uncertainties. 
Nevertheless, we see that even in the simple models considered here, the change of a in 
time can be compatible with both claims without any special fine-tuning. 

The meteoritic bound was discussed in detail in [26], where the bound —8 x 10~^ < 
Aa/a < 24 X 10~^ was derived. Note that in this case, the bound relates to the time 
average of Aa/a which is the quantity given in Table 1. Thus values in the table should 
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be restricted between -0.8 and 2.4. All models shown satisfy this constraint. 

The combined result of several recent atomic clock measurements limit the present day 
rate of change in the fine-structure constant to aja = (—0.9 ±4.2) x 10^^^ [45]. All models 
considered easily satisfy this bound. Finally, we compute the differential acceleration of the 
Earth-Moon system towards the Sun. The limit on Ag/g is 0.9 x 10"^^ [46]. As one can 
see, all of the models considered easily satisfy this bound as well. 

4.2 Early time evolution 

The experimental searches for variations in a provide a very useful constraint (or determi- 
nation, according to Ref. [21]) on the size of the quadratic term, ^ = An/(1 + b), where 
we dropped the subscripts in np and bp. Since ^ is typically 10~^ or smaller, one cannot 
determine whether this the result should be interpreted as 6 ^ 1 or n <^ 1 or both. 

However, there are two additional important observables, Aa at the cosmic microwave 
background epoch [47], and at BBN [42, 44, 48] that will provide more information on b 
and n. In fact, since the evolution of a at early times is monotonic, and since BBN and 
CMB data provide a similar level of sensitivity of order a few per cent to Aa/a, we will 
just use the BBN constraint as it extends much farther back in time. 

For both potentials, (3.4) and (3.5), we calculate the corresponding values of the field 
at the BBN redshift, ^ ~ 10^°. 

(j)BBN^-5.52 for \/(0) = \/oexp(A0V2) (4.4) 

(Pbbn ^ -11.9 for V(0) = 1/ocosh(A0) (4.5) 

These values of (psBN are calculated for the choice of A = 5, and can be used to estimate 
the size of the correction to the fine structure constant, Aa/asBN- Alternatively, we can 
use the scaling of the V{(f)) with a in Eq. (2.15) to get the same result. We conservatively 
assume that this quantity is limited at the level of 0.06 (see e.g. for the latest re-analysis of 
this bound [49]) based on the new derived ^He abundance of 1^ = 0.249 ± 0.009 [50]. The 
connection between V and Bp that we impose (3.2) does not allow for a linearized form of 
Bp{(f)) (3.1) at the high redshift associated with BBN, because the potentials we consider 
in this paper are very steep. 

Using Eq. (4.1) together with the ansatz (3.2), the constraint Aa/a{(f)BBN) becomes 

/ 1 -I- /) \ " 

' ^ >0.94 (4.6) 



b + Vi(l>BBN)/Vo 



This constraint is shown in Fig. 6 for our two choices of potentials, both with A = 5. In the 
light shaded region, both potentials satisfy the BBN bound on a. In the medium shaded 
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Figure 6: Constraints on the parametrized relation between B and V coming from the early- 
and late-time evolution of the scalar field. The constraints from the early-time evolution 
determine the allowed region in the b — n plane for the both potentials (light shaded) and 
for the cosh((^) potential alone (medium shaded). Both potentials fail the constraint in the 
dark shaded region. The early-time evolution constraint is shown by the near horizontal 
lines at low b. The values for ^ are taken from the table: ^ x 10^ = 5.8, 22, 310, and 1900 
are shown by the dashed, solid, dot-dashed, and dotted lines respectively. 

region, only the cosh(0) potential satisfies the bound. In the dark shaded region, both 
potentials lead to variations in a at the time of BBN which are in excess of the constraint. 
Also shown in Fig. 6 is the constraint obtained in the previous section from the late-time 
evolution of 0. Here we show the relation ^ = An/(1 -|- b) for the four values of ^ given in 
Table 1. The values of ^ x 10^ = 5.8, 22, 310, and 1900 are shown by the dashed, solid, dot- 
dashed, and dotted lines respectively. As one can see, when n < 10~^, the BBN constraint 
effectively disappears. At the same time, the choice of .^ ~ 10~^ and n > 10^'^ is clearly 
forbidden by the BBN constraint. Given this constraint on n, we can see that only very 
shallow exponential profiles of Bf{(J)) are allowed. 



5 Coupled quintessence 

In this section, we investigate changes to the cosmological evolution of the scalar field 
that can be induced by a non-trivial function Bm{4>)- The case of the so-called "coupled 
quintessence" has been addressed in the literature already [38], but was specialized to the 
particular case of an exponential potential with a linear coupling to matter. 

The function Bm{4>) expresses the dependence of the mass of any non-relativistic 
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particle, including dark matter particles. It is clear that if this mass experiences significant 
changes during the evolution of the Universe, the scaling of the matter energy density will 
differ from a~^, while the number density of dark matter particles n will continue to scale 
as a~^. The cosmological evolution of also changes, because is driven by an effective 
potential that now consists of two parts, 

\/(0)eff = V{(t)) + p™(0) = \/(0) + 5„(0)n. (5.1) 

It is not possible to say a prion which term is more important for the cosmological evolution 
of0. 

If Bm{<P) coincides with Bf{4>), we can use the insight gained from the previous section 
to conclude that Bm{4>) would provide no significant change on the cosmology from z = Oto 
zcMB- This is because the Aa/a constraints can be re-interpreted as a very tight bounds 
on np and hp-, and as a consequence on n™ and hm- However, Bm{(t>) and Bpiip) need 
not be the same. In fact, in simplest models (for example in the Brans-Dicke model, or 
string-inspired models with a modulus field 0), these functions typically are different. In 
the rest of this section, we assume that Bm{4>) and Bp{(l)) and not related and investigate 
possible constraints on rim and hm- 

For definiteness we choose hm = 1, and then Bm(4>) = (V^(0))"'", according to our 
parametrization. If we choose Um > 0, the effect of i?m(0) causes pm scale faster than a~^. In 
fact, if the dynamics of the scalar field is driven mostly by V{(j)) in the radiation dominated 
epoch, we can derive a simple upper bound on Um- Indeed, according to (2.15), the scaling 
of V^(0) in the radiation-dominated regime is a~*. Then it is clear that Bm{<P) ~ 1^(0)^^*^ 
would lead to Bm scaling as a~^, which also means that pm ~ a~'^. Thus, Um < 1/4, since 
for Hm > 1/4 the dark matter energy density would be redshifted in the same way or even 
faster than radiation. This would for example, greatly increase the energy density of dark 
matter during BBN. Therefore, Um < 0.25. In fact if we argue that the energy density in 
matter should not contribute more than the energy density in radiation at BBN, we can 
obtain a slightly stronger bound. For simplicity, let us take the ratio of matter to radiation, 
r, today to be 10*^ (a more accurate number would be ~ 5500). Then at the time of BBN, 
the ratio would be 

tbbn ~ roa^BBN (5-2) 

For rsBN ^ 1) ''"o ~ 10"^, and assN ~ 10^^°, we have n < 6/40. 

In Fig. 7 we show the changes in the cosmological evolution of the scalar field, matter and 
radiation that result from choosing Um = 0.1. As one can see, there are several important 
differences from the uncoupled case. Because of the change in the mass of the dark matter 
(as described above), matter domination occurs somewhat earlier. In this case, matter 
domination starts at a; ~ —9.6 or at z ~ 15000 and ends at a; ~ —0.33 corresponding to 
z ^ 0.4. Perhaps the most dramatic difference is seen in the evolution of the EOS parameter 
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Figure 7: a) The evolution of different components of the energy density as a function of x 
for the case of n^ = 0.1 and 6^ = 1? b) The evolution of under the same conditions. 

uj^. Instead of falling from the tracking value of 1/3 towards -1, as in the uncoupled case, 
here we see that cu^ increases towards +1, before ultimately dropping to —1 at later times. 
This is due to the fact that when the matter density increases, the effect on the potential 
V{(j))eg given in Eq. (5.1) drives the evolution of faster. The field begins to quickly roll 
and the energy density in becomes dominated by kinetic rather than potential energy. 
The evolution of is shown in Fig. 7b. As one can see, around x ~ —10, begins to 
change rapidly. Hence o;,^ ~ 1. (This is also the reason for the curious bump in Q^ seen in 
Fig. 7a.) However, ultimately approaches and settles at the minimum of the potential. 
At that time, u;,^ quickly moves towards -1 and the Universe becomes dominated by the 
dark energy. 



Finally we comment on the effect of the coupling to matter on the variation of a. The 
evolution of Aa is shown in Fig. 8 using the potential (3.4) with A = 5. One should note 
that because of the matter couphng, the field has moved quickly and at the present time 
has already passed the minimum twice for our choice of parameters. At late times {z < 3), 
the field is significantly closer to the minimum relative to the uncoupled case. As a result, 
the normalization using the QSO measurements results in significantly larger couplings, ( 
and/or ^. As one can see from the figure and Table 2, despite the fact that the field is 
closer to the minimum, the constraints from Oklo and the meteoritic abundances are in 
fact more severe. Only the results based on the normalization using Aa/a = 0.06 x 10~^ 
are compatible with these constraints. Similarly the differential acceleration is predicted 
to be much larger in this case, but the equivalence principle constraints are still satisfied. 
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In contrast, because the field has nearly settled to its minimum, and hence a is much 
smaller than in the uncoupled case. 

One can also determine the bounds on np and hp for this model. The late-time evolution 
constraint is very similar to, but slightly weaker than the result shown in Fig. 6 for the 
exp(A0^/2) case. The early-time evolution constraint is weaker than all of the cases shown 
in Fig. 6 and does not exclude any region of the h — n plane which is not already excluded 
by the late-time constraint. 



Table 2: As in Table 1 for the coupled case with B^ ^ 0, and bm = 1 and rim = 0.1. 



V(0) 


(^)3 


[ Aa\ 


e 


C 


[ Aa\ 

{ a J0.14 


I a JO.45 


f (yr-^) 




exp(A0V2) 


-5.4 


-5.7 





-150 


-0.85 


-1.3 


0.38 


-1000 




-5.4 


-7.4 


40000 





2.5 


2.6 


-0.034 


-13000 




-0.56 


-0.6 





-16 


-0.089 


-0.14 


0.039 


-11 




-0.44 


-0.6 


3200 





0.20 


0.21 


-0.0027 


-88 



6 Conclusions 

We have analyzed the cosmological evolution of the scalar field driven by its self-interaction 
potential, V{(f)), and its possible couplings to matter, Bm{<P)- Following an original Damour- 
Polyakov hypothesis, we assumed a common extremum for the functions V{(f)) and B„i{(f)) 
so that the evolution of the scalar field is driven towards that extremum, chosen here to 
be at = 0. With our choices of potentials, the evolution of the scalar field occurred 
in the tracking regime throughout the radiation dominated epoch. For definiteness we 
chose two potentials that ensure tracking and have a minimum, V{(f)) = Voexp(A0^/2) and 
V{(j)) = cosh(A0). 

We have seen that the coupling of the scalar field to the electromagnetic field is severely 
restricted by the late-time evolution of the field due to limits obtained from spectral mea- 
surements of quasar absorption systems, the Oklo phenomenon and meteoritic data. Models 
with quadratic couplings of the scalar field to radiation are typically more sensitive to the 
quasar absorption spectra- derived constraints than to the Oklo bound. The late-time evo- 
lution of the scalar field can be very nonlinear, and Aa/a can naturally be close to zero at 
certain redshifts, while having significant deviations over a range of the redshifts. Thus the 
experimental results, Refs. [21] and [22], claiming a non-zero result over a range of redshifts 
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Figure 8: As in Fig. 4 for the potential (3.4) with A = 5. In this case, Bm 7^ 0, and "we 
have chosen rim = 0.1 and 6m = 1- 
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0.5 — 3 and zero result at z = 1.5, do not have to be viewed as contradictory. We also find 
that the early-time evolution of the field severely constrains the form of the coupling -Bf(0) 
with respect to its relation to V {</)), due to BBN constraints on Aa/a. The constraints 
coming from BBN are complementary to those from the late-time change in a. 

Finally we have considered the case where the quintessence field is coupled to matter. 
Once again, BBN restricts the relationship between B and V due to the altered scaling of the 
matter energy density with the scale factor. Unlike the uncoupled case, the Universe briefly 
enters a period where the equation of state for quintessence is very stiff, before it evolves 
to a more standard dark energy EOS with u^ ~ —1. In order to have a O(10~^ — 10~®) 
relative change of a at redshifts 2; ~ 1, these models typically require larger values of 
coupling between the scalar field and electromagnetic radiation. 
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